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1 Introduction 

Integrable models of classical statistical mechanics [1, 2] and of two-dimensional quantum field theories 
[3], have a central object: the 7^-matrix TZ{u), where u is a spectral parameter, acting on the tensor 
product V (E)V oi a. given vector space V and being a solution of the Yang-Baxter (YB) equation 

ni2{u)ni3{u + v)n23{v) - 7^23(v)7^l3(w + v)tZi2{u), (i.i) 

in T/i «) T/2 ^ t/3, where 71x2 ^Tl(g) 1, 7^23 = f 7^, etc. 

The solution TZ{u) of Eq.(l.l) is said to be semiclassical if TZ{u) also depends on an additional pa- 
rameter T] in such a way that 

TZ{u,r]) = l + T]r{u)+o{T]^), (1.2) 

where 1 is the identity operator on the space V ^V. The "classical r-matrix" obeys the equation 

[^12 (m), 7-13 (u + v)+ r23{v)] + [ri3(M + v),r23{v)] = 0. (1.3) 

This equation, called the classical Yang-Baxter equation, plays an important role in the theory of classical 
completely integrable systems [4]. 

Nondegenerate solutions of (1.3) in the tensor product of two copies of a simple Lie algebra g , 
rij{u) e gi iX) gj , i,j ~ 1, 2, 3, were classified by Belavin and Drinfeld [5]. 

The classical YB equation has an interplay with conformal field theory, shortly described in the 
following way: in the skew-symmetric case rji{—u) + rij{u) = 0, it is the compatibility condition for the 
system of linear differential equations 

K =2^ry(z, -Zj)*(zi,...,ZAr) (1.4) 

in N complex variables zi,...,zn for vector- valued functions ^P with values in the tensor space V = 
V^ (E) ■ ■ ■ <E) V^ . K is a coupling constant. 

In the rational case [5], very simple skew-symmetric solutions are known: r{u) — C2/U, where C2 S 
g (g) g is a symmetric invariant tensor of a finite dimensional Lie algebra g acting on a representation space 
V. The above system of linear differential equations (1.4) is known as the Knizhnik-Zamolodchikov (KZ) 
system of equations for the conformal blocks of the Wess-Zumino-Novikov-Witten (WZNW) models of 
the conformal field theory on the sphere [6]. 

The algebraic Bethe Ansatz [7] was formulated in a parallel reasoning to the representation theory 
of the classical Lie algebras and it is a powerful method in the analysis of integrable models. Besides 
describing the spectra of quantum integrable systems, the algebraic Bethe Ansatz is also used to construct 
exact and manageable expressions for correlation functions [2] . Various representations of correlators were 
found by Korepin [8], using this method. 



The Babujian and Flume work [9] unveils a link between the algebraic Bethe Ansatz for the theory 
of the Gaudin models[12] and the conformal field theory of WZWN models. In their approach, the wave 
vectors of the Bethc Ansatz equation for an inhomogcneous lattice model render, in the scmiclassical 
limit, solutions of the KZ equation for the case of simple Lie algebras. For instance, in the su(2) example, 
the algebraic quantum inverse scattering method [7] allows one to write the following equation 

r(w|z)$(ui.....Up) = A(u, ui,-- • ,Up|z)$(ui,- •• ,Up) - V" — ^ . (1.5) 

Q — 1 

Here t(u\z) denotes the transfer matrix of the rational vertex model in an inhomogcneous lattice acting on 
an iV-fold tensor product of su(2) representation spaces. $" meaning $" — <&(mi, • • • Mq,_i,u, u„+i, ...,Up) 
; Ta{u\,--- ,Up|z) and A(u,'Ui,--- ,Mp|z) being c number functions. The vanishing of the so-called 
unwanted terms, Ta = 0, is enforced in the usual procedure of the algebraic Bethe Ansatz by finding the 
parameters ui, ..., u^. In this case the wave vector <I>(ui, • • • , Up) becomes an eigenvector of the transfer 
matrix with eigenvalue A(u,'Ui, • • • ,Up|z). If we keep all unwanted terms, i.e. Ta ^ 0, then the wave 
vector $ in general satisfies the equation (1.5), named in [10] as off-shell Bethe Ansatz equation (OSBAE). 
There is a neat relationship between the wave vector satisfying the OSBAE (1.5) and the vector- valued 
solutions of the KZ equation (1.4): The general vector valued solution of the KZ equation for an arbitrary 
simple Lie algebra was found by Schechtman and Varchenko [11]. It can be represented as a multiple 
contour integral 



^(zi, . . . , Zf^) — (p ■ ■ ■ d) X{ui, ..., Up\z)(j){ui, ..., Up\z)dui ■ ■ ■ dup. (1-6) 

The complex variables zi, ...,zi\j of (1.6) are related with the disorder parameters of the OSBAE . The 
vector valued function (f>(ui, ...,Up\z) is the scmiclassical limit of the wave vector $(ui, ..., Up|z). In 
fact, it is the Bethe wave vector for Gaudin magnets [12], but "off -shell ". The Bethe Ansatz for 
the Gaudin model was derived for any simple Lie algebra by Reshetikhin and Varchenko [13]. The 
scalar function X{ui, ...,Up\z) is constructed from the scmiclassical limit of the A(u = Zk]Ui, ...,Up\z) and 
the J-a(ui, ■ ■ ■ ,Up\z) functions. This representation of the iV-point correlation function shows a deep 
connection between the inhomogcneous vertex models and the WZNW theory. 

In this work we show that this idea also holds for the case of the semiclassical limit that corresponds 
to the sl{2\l)^'^' trigonometric r-matrix [15, 16]. 

The paper is organized as follows. In Section 2 we present the algebraic Bethe Ansatz for the sZ(2|l)'^^) 
vertex model. Here the inhomogcneous Bethe Ansatz is read from the homogeneous case previously 
derived for the 19- vertex models [14]. We also derive the off-shell Bethe Ansatz equation for this vertex 
model. In Section 3 , taking into account the scmiclassical limit of the results presented in the Section 
2, we describe the algebraic structure of the corresponding Gaudin model. In Section 4, data of the 



off-shell Gaudin equation are used to construct solutions of the trigonometric KZ equation. Conclusions 
are reserved for Section 5. 



2 Inhomogeneous Algebraic Bethe Ansatz 

Consider ^ = Vb ® Fl a Z2-graded vector space where and 1 denote the even and odd parts respectively. 
The components of a linear operator A (E) B in the graded tensor product space V CS) V result in matrix 
elements of the form 



(A^ByJ^ = (_)p(/3)(p(a)+p(7)) A^^Bf,s 

s s 

and the action of the permutation operator V on the vector \a) (^ \f3) G V (^ V is given by 

V \a) ® 1/3) = (-)P(")P(« 1/3) |a) =^ {V)l'^ = (-)p(")p('3)<5„^ 5,s. 

where p{a) ~ 1 (0) if \a) is an odd (even) element. 

Besides 7?, , it is usual to consider the matrix R — VTZ which satisfy 

Rl2{u)R2z{u + v)Ri2{v) ^ R23{v)Rl2{u + v)R23(u). 



(2.1) 



(2.2) 



(2.3) 



The regular solution of the graded YB equation of the 19-vertex model sZ(2|l)'^^) is given by [15, 16] 
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Consider now an inhomogcncous vertex model, where to each vertex we associate two parameters: 
a global spectral parameter u and a disorder parameter z. In this case, the vertex weight matrix TZ 
depends on u — z and consequently the monodromy matrix defined bellow will be a function of the 
disorder parameters Zi. 

The graded quantum inverse scattering method is characterized by the monodromy matrix T(u\z) 
satisfying the equation 



R{u-v) T{u\z)(g)T{v\z) = T{v\z)(g)T{u\z) R{u - v) 



(2.6) 



whose consistency is guaranteed by the graded version of the YB equation (1.1). T{u\z) is a matrix in 
the space V (usually called auxiliary space) whose matrix elements are operators on the states of the 
quantum system (which will also be in our work the space V). The monodromy operator T{u\z) is defined 
as an ordered product of local operators £„ (Lax operator), on all sites of the lattice: 



T{u\z) = Cn{u- zn)Cn-i{u- zn-i) ■ ■■ Ci{u~ zi). 



(2.7) 



We normalize the Lax operator on the n quantum space to be given by 
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(2.8) 



Ll(u)^ a, (3 ~ 1, 2, 3 are 3 by 3 matrices acting on the n site of the lattice, so the monodromy matrix 
has the form 



T{u\z) 



where 



T,j{u\z) 



hJ 



1,2,3. 
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(2.10) 



The vector in the quantum space of the monodromy matrix T{u\z) that is annihilated by the operators 
Tij(u\z), i > j {C'i{u\z) operators, i = 1,2,3) and it is also an eigenvector for the operators Tii(u\z) ( 
Ai(M|z) operators, i = 1, 2, 3) is called a highest weight vector of the monodromy matrix T{u\z). 

The transfer matrix t(u\z) of the corresponding integrable spin model is given by the supertrace of 
the monodromy matrix in the space V 

3 

t{u\z) = ^(-ir« Tu(u\z) = Ai{u\z) - A2{u\z) + A3(u\z). (2.11) 

The algebraic Bethe Ansatz solution for the inhomogeneous sl{2\l)'-^^ vertex model can be obtained 
from the homogeneous case [14]. The proper modification is a local shift of the spectral parameter 

U ^ U — Zi. 

Here we will define the functions used in our algebraic Bethe Ansatz: 



Xi{u) sinh(u + 277) 2:3 (w) cosh(?i — 77) 

z[u) = — i—r = r-r , ?/(« - - 



uj{u) 



X2{u) sinhw ' xq{u) sinh2r7 

xi{u)x3{u) sinh(u + 277) cosh(u — 77) 

X4{u)x3{u) — xe{u)xQ{u) sinh(7i — 27;) cosh(u + 77) ' 




Ziu,-u,) = I ^("'=-"^) ^ ''I ^>{ . (2.12) 

[_ Z[Uk — Uj)Ll)[Uj — Uk) II K < J 

We start defining the highest weight vector of the monodromy matrix T(u\z) in a lattice of N sites as 
the even (bosonic) completely unoccupied state 

/ 1 \ 

(2.13) 

Using (2.10) we can compute the normalized action of the monodromy matrix entries on this state 
A{u\z)\0) = X,{u\z)\0), C,(w|z) |0) = 0, B,iu\z)\0)^{0,m, 

XMz) = li ( V ^=1-2,3. (2.14) 

^^ X2{U- Za) 

The Bethe vectors are defined as normal ordered states ^„(ui, • • • , 7i„) which can be written with aid 
of a recurrence formula [17]: 

^niui,...,Un\z) = Si(7ii|z)^„_i(M2,...,M„|z) 

-B2iu,\z)y] -^l("j-|^) TT Z{uk-U,)^n-2{u2,...,U,,...,Un\z), (2.15) 



with the initial condition ^Pq = |0) , ^i(mi|z) = Bi{ui\z) |0). Here Uj denotes that the rapidity Uj is 
absent: "^{ujlz) ~ "^{ui, ...,Uj_i,Uj^i, ■ ■ ■ ,u„|z). 

The action of the transfer matrix t(u\z) on these Bethe vectors gives us the following off-shell Bethe 
Ansatz equation for the s/(2|l)*^^-' vertex model 



n n j — 1 

j=i i=2 1=1 (2.16) 

We now briefly describe each term which appear in the right hand side of (2.16) ( for more details 
the reader can see [14]): in the first term the Bethe vectors (2.15) are multiplied by c-numbers functions 
A„ = A„(u,wi,...,u„|z) given by 

k=i ^J-^u;{u-Uk) ^^\xs{u-uk) (2.17) 

The second term is a sum of new vectors 

multiplied by scalar functions J^^''^ given by 

n n 

^j""'' = Xiiuj\z) Y[ Z{uk - Uj) + {-TX2{u.j\z) J] Z{u.j - Uk). (2.19) 

Finally, the last term is a coupled sum of a third type of vector-valued functions 

*li-2) = B2iu\z)^'n^2{ui,Uj), (2.20) 

with coefficients 

n 
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where Gij , Y/j , Fij and Hij are additional functions defined by 
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In the usual Bethe Ansatz method, the next step consist in impose the vanishing of the so-called unwanted 
terms of (2.16) in order to get an eigenvalue problem for the transfer matrix. 

We impose ^ — and J^j^" — into (2.16) to recover the eigenvalue problem. This means that 

^„('Ui, ...,Un\z) is an eigenstate of r(u|z) with eigenvalue A„ , provided the rapidities Uj are solutions of 
the inhomogeneous Bethe Ansatz equations 

J = l,2,...,n. (2.23) 

3 Structure of the sU2|l)^^^ Gaudin Model 



In this section wc will consider the theory of the Gaudin models. To do this we need to calculate the 
semiclassical limit of the results presented in the previous section. 

First, we recall that the superalgebra sZ(2|l) ~ sZ(l|2) is the {N — 2) extended supersymmetric version 
of sl{2) and contains four even (bosonic) generators iJ, W,X^ which form the Lie algebra sl{2) © U{1) 
and four odd (fermionic) generators V^ , W whose non- vanishing commutation relations read as [18, 19]: 



[H,X±] = 


= ±^^ 


[X+,X-] ^2H, 


{v+,v-}- 


--i" 


[H.V^] - 


- 4-^ 


[X^,V^] ^v^, 


{t/±,y±} = 


--¥ 


[x^.v^] - 


= 0, 


[X^,W^] =0, 


{V^,W^}-^ 


= 0, 


[A±,iyT] = 


= ^^ 


[w,w^]^tIx^, 


[W,X^] = 


4"- 


{V^,W^} -- 


= ±^^^ 


{w+,w-} = Ih, 


{W^,W^} = 


= A±. 



(3.1) 






/ 1 


= 







\ . 


l/ 


'000 


-10 


2 
\ 


1 



The quadratic Casimir operator is 

C2=H'^ -W"^ + X-X+ + 2W-W"^ -2V-V"^, (3.2) 

and the elementary representation is three-dimensional, given by 

. / 1 \ , / 1 

W = - 2 , H=-[ 

"^Vooi/ ^\oo-i 

X- = I I , y+ - i 1 I , V- 

"^ \ 

,/010\ ,/000 
W^+ = - -1 , VK" = - 1 I (3.3) 

^\000/ ^\010 

Here the basis is {|1) , |0) , | — 1)}, where the first and third vectors are considered as even and the second 
as odd, i.e., the grading is BFB. 

In order to expand the matrix elements of T{u\z), up to an appropriate order in 77, we will start by 
expanding the Lax operator entries defined in (2.8): 

j(n) ^ ^ 2->V„ + 7^„cosh(2u-2^„) 2/^1 2 1 nl-Un 

" "^ sinh(2u-2z„) "^ ^2 " 4cosh(u-z„)2 
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^(„) ^ 1 + 2,^ - ^" .- ^" -y- - ^-") + V f ^^2 - 1 ^j + ^" ,, I + o(,-). 
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(3.4) 



and for the elements out of the diagonal we have 



12 ' sinh(2u-2z„) ^ ^'^' 21 / sinh(2u-2z„) 
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where V± = 2t/±, A:'± == 2X± , >V± = 2VK±, H = 2iJ and W = 2VK. 

Substituting (3.4) and (3.5) into (2.10) we will get the semiclassical expansion for the monodromy 
matrix entries. For the diagonal entries we get 

yli(7i|z) = l + 277Af^(7i|z)+4772Af)(7i|z) + o(7;3), z = 1,2,3. (3.6) 



where, the first order terms are given by 
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For ofF-diagonal elements we only need to expand them up to the first order in rj 
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From these we have the following expansion for the transfer matrix (2.11): 
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Now we will consider the second order term of t{u\z): 
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which, with aid of the identity 
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sinh(2M — Iza) sinh(2u — 2z(,) sinh(22;o — 2zh) sinh(2M — 2za) sinh(2u — 2zb) ' 
can be written in the form 



TV 



-'''H-)-Y.i^- 



Qa{u) 



N 



^^^ e2"-2za sinh(2u - 2zq) ^^ " 2 cosh(u - z^) 

a—l a—l 



E(«a 



1 1 



(3.12) 



(3.13) 
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where 

Gain) = J2-r-u^ — --—{2 + cosh{'2u-2za)cosh{2u-2zb)Hai)nb 

-Wa ®Wb + ^ [x+ I. X^- + X^ A'+) 

+ 0--"+^^ (W- ^ W+ + V+ ® V,-) - e-»--' (W+ ® W- + V- i> V+)} } 



(3.14) 



Here we observe that Gaiu) is nothing but the sum of semiclassical trigonometric r- matrices. This fact 
follows from the construction of the semiclassical r-matrices out of the quadratic Casimir: 

C2 = H^ -W^ + -{X'X+ + X+X-) + {W-W+ + V+V-) - [w+w- + V-V+). 

^ (3.15) 

The Gaudin Hamiltonians are defined as the residue of t{u\z) at the point u = Za- This results in N 
non-local Hamiltonians 

N 

G. = E-h(^I-^{2 + cosh(2.„-2..)7^.^H. 

-Wa cl Wb + i (^X+ ® X^ + X- ® A'+) 

+ c-^^+^^W- ® W+ + V+ cl V,-) - c^"-^HW+ ® W,- + V- ® V+)} , 
a = 1,2, ...,N. (3.16) 

satisfying 

Y^Ga^O, ^ = ^, [G,,Gfc] = 0, Wa,b. (3.17) 

a— 1 

In the next section we will use the data of the Bethe Ansatz presented in the previous section in order 
to find the exact spectrum and eigenvectors for each of these A^ — 1 independent Hamiltonians. 

We complete this section deriving the sl{2\l)'-^^ Gaudin algebra from the semiclassical limit of the 
fundamental commutation relation (2.6): The semiclassical expansions of T and R can be written in the 
following form 

T{u\z) = 1 + 2tj[1{u\z) + (3(u\z)] + o{rf), R{u) = V[l^ 2rj[r{u) + P{u)] + o{rf)\ . 

(3.18) 

where 

N 2 

/^H-)-E 3inh(2.-2..) - (^-1^) 
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Using (3.8-3.9) one can see that the "classical Z-operator " has the form 

/ W{u\z) + n{u\z) W- {u\z) - V- {u\z) X-{u\z) \ 



l{u\z) = 



W+{u\z) + V+{u\z) 



2W{u\z) W-{u\z) + V-{u\z) 

-W+{u\z) + V+{u\z) W{u\z)-n{u\z) J 



where 



N 



N 



-Wa 



H{u\z) - ^coth(2u-2z,)Ha, M«k)-5I-i.^o _o ^ 

a— 1 ^ ^ 



a=l 

N 



N 



W-{u\z) = J2- 



^-^ sinh(2ii — 2za) 

a—1 ■ ^ 



W-, W+{u\z)=J2- 



^-^ sinh(2w — 2za 

a—l ^ 



-w; 



N 



V~{u\z) = J2- 



-U + Za 



N 



^-^ sinh(2M — 2za) 

a— 1 ^ ' 

N 



V-, V+{u\z)^J2- 



^-^ sinh('2w — 2za 

a—l ^ 



■v^ 



Z_— / Qir 



X: 



N 



X+(u\z)^J2- 



x+ 



^-^ sinh(2'u — 2za) 

a—l ^ ' 



^-^ sinh(2M — 2za) ' 

a—l ^ ' 

The corresponding semiclassical r-matrix has the form 

riu) = — \ — jcosh2MH(|7i:- W®W+-(A:'+(|a'-+A'- |)A'+) 
sinh2u [_ 2 

+e" (W" ® W+ + V+ I) V") - C-" (W- ®YJ"^ ^V- ® V+). 



(3.20) 



(3.21) 



(3.22) 



Here we notice that (3.22) is equivalent to the r-matrix constructed out of the quadratic Casimir (3.15) 
in a standard way [20]. 

Substituting (3.22) and (3.20) into (2.6), we have 



Vl[u\z) ® l{v\z) + Vr{u - v) l[u\z) (g) 1 + 1 (g) l{v\z) 
l(y\z)®l(u\z)V + l(y\z)®\'\^\®l{u\z) Pr{u-v), 



(3.23) 



whose consistence is guaranteed by the graded classical YB equation (1.3). 

From (3.23) we can derive (anti-)commutation relations between the matrix elements of l(u\z). This 
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gives us the defining relations of tlie s^(2|l)(^' Gaudin algebra 



[n{u\z),H{v 

[X'^{u\z),W{v 

[X^{u\z),V^{v 

[W^{u\z),W{v 

[V^{u\z),X^{v 

[W^{u\z),X^{v 

[V=^(u|z),W(w 

[X^{u\z),n{v 

[v^{u\z),n{v 

[W^{u\z),Hiv 

{V^{u\z),V^{v 

{W'^{u\z),W'^{v 

{V^{u\z),W'^{v 

{W^{u\z),W^iv 

{V=^(u|z),V=^(w 



1^) 



= 0, [W(u\z),H{v\z)]^0, [W{u\z),W{v\z)]^0, 

= 0, [X'^{u\z),X^{v\z)] = 0, [X^{u\z),W^{v\z)]^0, 

= 0, {V^{u\z),W^{v\z)}^0, 



= ±4 



n{u\z) - H{v\z) 



sinh(2u — 2v) 

V^iu\z) -e±(-"+")V±(w|;j) 

sinli(2'u — 2v) 

V=F(u|z)-e±(«-")V=F(w|z) 

sinh(2M — 2v) 
>V=f(m|z) - e±(-"+^)W=F(w|z) 
sinh(2u — 2v) 
>V±(m|z) - e±("-^)W±(w|z) 



= ±2 



± 



± 



sinh(2u — 2v) 
cosh(2M - 2u)A'±(u|z) - A:'±(w|z) 
sinh(2M — 2w) 
cosh(2M - 2v)V^{u\z) - e±("-^)V±(w|z) 

sinli(2'u — 2v) 
cosh(2u - 2v)W'^{u\z) - c±(-"+^)W±(i;|z) 



z)} = - 
z)} = ± 



sinh(2u — 2v) 
eMu^v) (T-Clulz) - niv\z)) 
sinh(2M — 2v) 
e±(-''+") (n{u\z) - niv\z)) 






± 



T 



sinli(2M — 2w) 
e±(-"+")>V(u|z) - e±("-^)>V(w|z) 

sinh(2u — 2w) 
e±("-^)A'±(w|z) - c±(-«+^)A'±(w|z) 

sinli(2'u — 2v) 
e±(-"+")A:'±(M|z) - e±("-^)A:'±(w|z) 



sinh(2u — 2v) 
A direct consequence of these relations is the comniutativity of t^'^'{u\z) 

from which the commutativity of the Gaudin Hamiltonians Ga follows immediately. 



(3.24) 



(3.25) 
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4 Off-shell Gaudin Equation 

In order to get the semiclassical limit of the OSBAE (2.16) we first consider the semiclassical expansions 
of the Bethe vectors defined in (2.15), (2.18) and (2.20): 

*„(ui,...,u„|z) = (277)"$„(wi,...,u„|z) + o(ry"+i), 

*^- - 2(2:,)"+^ [^'("'"^^'"^"^'^'("'"^^"'"^^ t .((l-|z) + o(,"+2) 

^(„-i) - A2r?) sinh(2« - 2^*,) >^n-i{u,\z) + o[7^ ), 

*|l-2) = (2r?)"-iA'-(^i|z)<I>„_2(^;,^,k) + o(r;"), (4.1) 



where 



$„(mi,...,w„|z) = [W"(ui|z) - V"(ui|z)] $„-i(m2,---,u«|z) 

-A'"(ui|z) V -^ -$„_2(M2,Uj,Un|2), (4.2) 

^ cosh(ui — Uj) 

with $0 == |0) and $i(ui|z) = [W-{ui\z) - V-{ui\z)] $o- 

The corresponding expansions of the c-number functions presented in the OSBAE arc: (2.16) 

A„ = l + 2r7Ali)+VAi2)+o(ry3), (4.3) 



^j"-i) = 2r;(-)^+i/j"-^)+o(ry2), (4.4) 

/ \l+1 ( J-(n-l) f("-l) ~1 

-(n-2) _ „.„ v3 (-) J /^ Jj I , ,,/„4n 



'^ ^ '^ cosh(u/ - Mj) I sinh(2M - 2u0 sinh(2u - 2uj) j ^''' 



where 



(4.5) 



^i"^ - E • .J o , (4-6) 

^ sinh(2?i — 2za) 
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2^-^ cosh(u - ZaY 

a— 1 ^ ' 



and 



N n 

— > y {coth(u — Za) coth(u — Uj) + tanh(?i — Za) tanh(?i — Mj)} 

a=l j=l 

N 

+ y {coth(u — Za) coth(u — Zf,) + tanh(u — z^) tanh(u — z^)} 

a<b 

n 

+ 2, {coth(u — Uj) coth(u — Uk) + tanh(u — Uj) tanh(u — Uk) 
j<k 

4 
sinh(2u — 2uj) sinh(2M — 2uk) 



N 



(4.7) 



fj = X! coth(uj - Za) - ^ tanh(uj - Uk). (4.8) 

a=l k^j 

Substituting these expressions into the (2.16) and comparing the coefficients of the terms 2(2?7)"+^ 
we get the first non-trivial consequence for the semiclassical limit of the OSBAE: 

n 2f*-"^"'"-'9-' 

r(2)(«|z) $„(«!,.. .,u„|z) = Ai2) ^^{u^,...,Un\z)~Y.(^y .' <"'^; . (4.9) 

■'^-^ smh 2u — 2m,- 

Note that in this limit the contributions from ^/„_]^-, and ^Zj_2') ^'"^ combined to give a new vector 
valued function 

e^„_i) = [W-{u\z)e-^+"^ -V-{u\z)e"'^^] <^,,^,{u,\z) 

-X-{u\z) V —^ -$„_2((lj,Ufc|z), (4.10) 

where fc = k + 1 for /c < j and k ~ k for fc > j. 

Finally, we take the residue of (4.9) at the point u = Za to get the off-shell Gaudin equation: 



Ga^n{ui,...,Un\z) = ga$„ (ui , . .. , U„ |z) + V'(-)' — 



2f("-i)j,/ 



^^^ sinh(2'Ui - 2za) ' 
a = 1,2, ...,7V (4.11) 



(2) 

where ga is the residue of A„ 



N n 

9a = reSu^z^A^^^ = ^ coth(za - ^b) - X! coth(zQ - ui), (4.12) 

fc^iia 1 = 1 
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and 4>\n-i) is the residue of ©L^d 



aJ 



res,. 



.e^' 






The equation (4.11) allows us solve one of the main problem of the Gaudin model, i.e., the determi- 
nation of the eigenvalues and eigenvectors of the commuting Hamiltonians Ga (3.15): ga is the eigenvalue 
of Ga with eigenfunction $„ provided ui are solutions of the following equations / =0, i.e.: 

n N 

y^ tanh(Mj — Mfc) — /^ coth(Mj — Zg), j — l,2,...,n. (4-14) 

k^j a=l 

Moreover, as we will see in the next section, the off-shell Gaudin equation (4.11) provides solutions for 
the differential equations known as KZ equations. 

5 Knizhnik-Zamolodchickov equation 

The KZ differential equation 

«^^ = G.(z)vI/(z), (5.1) 

appeared first as a holonomic system of differential equations of conformal blocks in a WZW models 
of conformal field theory. Here '^{z) is a function with values in the tensor product Vi (8) • • • ® Vat of 
representations of a simple Lie algebra, k = k + g , where k is the central charge of the associated 
Kac-Moody algebra, and g is the dual Coxeter number of the simple Lie algebra. 

One of the remarkable properties of the KZ system is that the coefficient functions Gi{z) commute 
and that the form uj = J^i Gi(z)dzi is closed [13]: 

f = f, lft,6y = 0. (.2, 

Indeed, it was indicated in [13] that the equations (5.2) are not just a flatness condition for the form to 
but that the KZ connection is actually a commutative family of connections. 

In this section we will identify Gi with our s^(2|l)"' Gaudin Hamiltonians Ga derived in the previous 
section 

N -. 
Ga = y^^-TTT^ -—{cOsh{2Za-2zb)Ha^nb 

f-f smh{2za - 2zb) ^ 

+ e-^"+-'(W- ® W+ + V+ ® V,-) - e^"--HW+ ® W," + V" ® V+)} , 
a = l,2,...,iV. (5.3) 
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and show that the corresponding differential equations (5.f ) can be solved via the off-shell Bethe Ansatz 
method. 

Let us now define the vector-valued function ^(zi, ..., zn) through multiple contour integrals of the 
Bethe vectors (4.2) 



*(zi, ..., zn) — ([>■■■([> X{u\z)^n{u\z)dui...dun, (5.4) 

where X {u\z) ~ X (ui, ..., u„, zi, ..., z^r) is a scalar function which in this stage is still undefined. 
We assume that ^(zi, ..., zm) is a solution of the equations 

K^5%^^^^ = G,vI/(zi,...,z^), a- 1,2, ...,7V (5.5) 

where Ga are the Gaudin Hamiltonians (3.15) and k is a constant. 
Substituting (5.4) into (5.5) we have 

^3*fe_^ . / {-.^*,.(»l«) + ^nu\.f-^] <lu. (5,6) 

where we are using a compact notation f {o} du = § . . . § {o} dui ■ ■ ■ dun- 

Using the Gaudin algebra (3.24) one can derive the following non-trivial identity 

dza j-(' ' dui \sxnh{2ui-2za))' ^ ■ ' 

which allows us write (5.6) in the form 

a* _ // dX{u\z) " dX{u\z) ( 20|„_^) 



n 



!,,,„,), 5;,_y.-^_g.H_U.„ 



dza J 1 dza ^ dui \ sinh(2u( - 2za) 

X{u\z) . ^t^'^-'l A du. (5.8) 

smh(2M; - 2Za) / 



4ri^^i: 



It is evident that the last term of (5.8) vanishes, because the contours are closed. Moreover, if the scalar 
function A'(m|z) satisfies the following differential equations 

we are recovering the off-shell Gaudin equation (4.11) from the first term in (5.8). 

Taking into account the definition of the scalar functions /-" (4.8) and ga (4.12), we can see that 
the consistency condition of the system (5.9) is insured by the zero curvature conditions 9/ /dza = 

dga/duj. Moreover, the solution of (5.9) is easily obtained 

N n N n 

X{u\z) = Yl sinh(z„ - Zb)^/^ ]J cosh(Mj - Uk)^/^ J| J| sinh(za - %)"^/'*. (5.10) 

a<b j<k a—lj—1 

This function determines the monodromy of ^(zi, ..., z^r) as solution of the trigonometric KZ equation 
(5.5) and these results are in agreement with the Schechtman-Varchcnko construction for multiple contour 
integral as solutions of the KZ equation in an arbitrary simple Lie algebra [11]. 
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6 Conclusion 

In this paper a graded 19-vertex model was used to generalize previous rational vertex models results 
connecting the Gaudin magnet models to the semiclassical off-shell Bcthe Ansatz of these vertex models. 

Using the semiclassical limit of the transfer matrix of the vertex model we derived the trigonometric 
sl{2\l)^'^' Gaudin Hamiltonians. The reduction of the off-shell Gaudin equation to an eigenvalue equation 
gives us the exact spectra and eigenvectors for these Gaudin magnets. Data of the off-shell Gaudin 
equation were used to show that a hypergeometric type integral (5.4) is solution of the trigonometric KZ 
differential equation. 

In fact, this method had already been used with success to constructing solutions of trigonometric KZ 
equations [21, 22, 23] and elliptic KZ-Bernard equations [24], for the six-vertex model and eight-vertex 
model, respectively. 
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